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Abstract 

The synchronization transitions in Newman-Watts small-world neuronal networks (SWNNs) induced by time delay t and 
long-range connection (LRC) probability P have been investigated by synchronization parameter and space-time plots. Four 
distinct parameter regions, that is, asynchronous region, transition region, synchronous region, and oscillatory region have 
been discovered at certain LRC probability P= 1.0 as time delay is increased. Interestingly, desynchronization is observed in 
oscillatory region. More importantly, we consider the spatiotemporal patterns obtained in delayed Newman-Watts SWNNs 
are the competition results between long-range drivings (LRDs) and neighboring interactions. In addition, for moderate 
time delay, the synchronization of neuronal network can be enhanced remarkably by increasing LRC probability. 
Furthermore, lag synchronization has been found between weak synchronization and complete synchronization as LRC 
probability P is a little less than 1.0. Finally, the two necessary conditions, moderate time delay and large numbers of LRCs, 
are exposed explicitly for synchronization in delayed Newman-Watts SWNNs. 



•0-PLOS I ONE 



Citation: Qian Y (2014) Time Delay and Long-Range Connection Induced Synchronization Transitions in Newman-Watts Small-World Neuronal Networks. PLoS 
ONE 9(5): e96415. doi:10.1371/journal.pone.0096415 

Editor: Matjaz Perc, University of Maribor, Slovenia 

Received March 3, 2014; Accepted April 7, 2014; Published May 8, 2014 

Copyright: © 2014 Yu Qian. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted 
use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was supported by the National Natural Science Foundation of China (Grant Nos. 11105003), and the Natural Science Foundation of the 
Education Bureau of Shaanxi Province of China (Grant No. 201 3JK061 9). The funders had no role in study design, data collection and analysis, decision to publish, 
or preparation of the manuscript. 

Competing Interests: The author has declared that no competing interests exist. 
* E-mail: qianyu0272@163.com 



Introduction 

Synchronization phenomena are common in nature and can be 
extensive observed in various realistic systems, especially in 
neuronal networks, biological systems and ecological systems [1,2]. 
Synchronization has been widely studied both theoretically and 
experimentally for decades. Several kinds of synchronization have 
been discovered in theoretical researches, such as complete 
synchronization, weak synchronization, lag synchronization, phase 
synchronization and generalized synchronization [3-7]. Complete 
synchronization indicates the coincidence of states of coupling 
systems, Xi(/) = X2(0 [3]. Lag synchronization described in 
Ref. [5] means the coincidence of shifted in time states of two 
systems, Xi(/ + tq) = X2(0- Experimental studies have shown that 
synchronous oscillations can emerge in many special areas of 
brain, especially in olfactory system or hippocampal region [8- 1 0] . 
In recent years, synchronization in neuronal networks and brain 
systems has attracted much attention. Synchronous oscillations in 
these systems are related to some specific and important 
physiological functions, such as olfaction [11], visual percep- 
tion [12], cognitive processes [13], and information processing [14]. 

Recently "small-world" network has been proposed by Watts 
and Strogatz, which takes into account both local and long-range 
interactions [1 5] . It is found that the existence of a small fraction of 
long-range connections (LRCs) can essentially change the features 
of the given systems [16-19]. These LRCs do exist in neuronal 
networks and do play crucial roles in deciding the specific 



physiological functions. The interactions from the long-range 
connected neurons must be time delayed due to the finite 
propagation velocities in the conduction of signals along neuron 
axons [20]. And the effects of time delays on self-organized 
spatiotemporal dynamics in neuronal systems have been exten- 
sively investigated. Lots of interesting phenomena have been 
discovered in recent decades [2 1-36] . For example, Dhamala et al. 
have investigated the enhancement of neural synchrony by time 
delay [21]. Ko et al have found that time delay can destabilize 
synchronous states and induce near-regular wave states [23]. 
Significantly, Wang et al. have discovered that time delays can 
enhance the coherence of spiral waves [27], tame desynchronized 
bursting [28], induce stochastic resonances [29] and synchroniza- 
tion transitions [30-32] , and can cause synchronous bursts [33] and 
complex synchronous behavior [34]. Moreover, Yu et al. have 
demonstrated the synchronization transitions in delayed neuronal 
networks with hybrid synapses [35,36]. Although remarkable 
advances have been achieved in the field of delayed neuronal 
networks, the underlying mechanisms behind time delay induced 
spatiotemporal dynamic and related synchronization transitions 
are far from being fully understood. In addition, the lag 
synchronization, to our knowledge, has not been identified in 
delayed neuronal network. These are the tasks we aim to explore. 

In this paper we extend the subject by systematically investi- 
gating time delay and long-range connection induced synchroni- 
zation transitions in Newman-Watts small-world neuronal net- 
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works (SWNNs). By introducing synchronization parameter and 
plotting spatiotemporal patterns, four distinct parameter regions, 
i.e., asynchronous region, transition region, synchronous region 
and oscillatory region, have been found at certain LRC probability 
P= 1.0. Interestingly, desynchronization and oscillating behaviour 
of the order parameter are observed in oscillatory region. More 
importantly, the mechanisms of synchronous oscillations and the 
transition from non-synchronization to complete synchronization 
are discussed. Moreover, we consider the spatiotemporal patterns 
obtained in delayed Newman-Watts SWNNs are the competition 
results between long-range drivings (LRDs) and neighboring 
interactions. A new order parameter, LRD proportion, is used 
to verify our point of view. And the four distinct parameter regions 
can also be revealed by LRD proportion clearly. In addition, for 
moderate time delay, the synchronization of neuronal network can 
be enhanced remarkably by increasing LRC probability. Further- 
more, lag synchronization has been found between weak 
synchronization and complete synchronization as LRC probability 
P is a little less than 1.0. And the mechanism is revealed. Finally, 
the two necessary conditions, moderate time delay and large 
numbers of LRCs, are exposed explicitly for synchronization in 
delayed Newman-Watts SWNNs. 

Mathematical Model and Setup 

We start from a one-dimensional (ID) regular ring that 
comprises TV =100 identical excitable Bar-Eiswirth neurons [3 7] 
with periodic boundary condition, and each neuron has two 
nearest neighbors. The evolvement of the ID neuronal network is 
governed by the following equations: 



dUi{t) 1 V/(0 + ^ , 

= - -Ui{t)[Ui{t)-\\[Ui{t) — ] 

+ D[ui-x{t) + Ui^x{t)-2ui{t% 



dvijt) 
dt 



=f[Ui{t)\-Vi{t\ 



(1) 



(2) 



where /= 1,2, . . . ,7V. The function f\ui{t)\ takes the form: 
f[Ui{t)\=0 for t/.(0<J; f[Ui{t)\ = \-6.15Ui{t)[Ui{t)-\f for 

\^<Ui{t)<\\ and f[Ui{t)\ = \ for Ui{t)>\. Here variables u and 

V are the activator and inhibitor variables, respectively. The small 
relaxation parameter e represents the time ratio between activator 
u and inhibitor v. The dimensionless parameters a and b denote 
the activator kinetics with b effectively controlling the excitation 
threshold. D is the coupling intensity which decides the interaction 
strength between neighboring neurons. The system parameters are 
kept throughout this paper as a = 0.84, /? = 0.07, e = 0.04 and 
i) = 0.5. Therefore, the local dynamics can describe typical 
excitability of neurons where u represents the membrane potential, 

V is the somatic inhibitory current. The diffusive couplings 
simulates electrical conjunction interaction between neurons. 

Based on the ID periodic regular ring, we construct delayed 
Newman- Watts SWNNs [38] by introducing LRCs such that each 
neuron receives an unidirectional time delayed LRD from a 
randomly chosen cell with probability P [39,40] . We thus add an 
additional coupling term to Eq. (1) if neuron / receives an 
unidirectionally time delayed LRD from cell j. Now the delayed 
Newman- Watts SWNN is governed by the following equations: 



dui{t) 
IF 



= - - Ui(t)[Ui{t) - 1] [Ui(t) - ] 

+ D[Ui-x{t) + Ui^x(t)-2Ui{t)\ 

+ D[uj{t-x)-Ui{t% 



(3) 



dt 



=f[um-yiity 



(4) 



In Eq. (3) cell j is randomly chosen in the ID periodic regular 
ring and T is the time delay in information transmission. By 
manipulating LRC probability P, we can obtain different kinds of 
time delayed Newman-Watts SWNNs. The schematic diagram of 
the considered networks for different LRC probability with 10 
neurons is illustrated in Fig. 1 . Here we should mention that for a 
given LRC probability there are a lot of network realizations. For 
a specific network structure, the interactions between neighboring 
neurons are bidirectional (shown by bidirectional arrowed lines), 
while the LRDs are unidirectional (shown by unidirectional 
arrowed lines). Time delays are only considered in these 
unidirectional LRDs, which will cause inhomogeneity in informa- 
tion transmission between neighboring and long-range interac- 
tions. As we know that the interactions from neighboring neurons 
are usually instantaneous in actual biological systems. And the 
LRDs from distant cells will have time delays due to the finite 
propagation velocities. Therefore, the model considered in present 
paper may be more realistic, and the results obtained may be more 
practical. 

In this paper, the delayed Newman-Watts SWNNs are 
integrated by forward Euler integration scheme with time step 
A/ = 0.001. The initial variables (w/(/ = 0),V/(/ = 0)) are randomly 
given between 0 and 1 for each simulation. To investigate the 
synchronization transitions in delayed Newman-Watts SWNNs 
quantitatively, the synchronization parameter R will be used, 
which has been introduced in the previous study [41]. It is 
numerically calculated as: 



R- 



where 



<u{tf>- <u(t)>^ 



1 ^ 



(5) 



(6) 



The angular brackets denote the average over time. In present 
paper the synchronization parameters are calculated over last 30 
time units. From Eq. (5) it is evident that the larger the 
synchronization parameter R is, the more synchronization is 
realized in neuronal network. Accordingly, the value of R close to 
unity indicates all neurons in the network are in complete 
synchronization. Therefore, the synchronization parameter R is 
an excellent indicator to reveal the spatiotemporal synchronization 
in delayed Newman- Watts SWNNs and the related transitions. 

To guarantee the statistical accuracy with respect to the network 
structure and initial condition, 10 independent samples are 
executed for each set of parameter values in the simulation. And 
we will use 
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Figure 1. Schematic diagram of the considered Newman-Watts small-world neuronal networks for different long-range connection 
(LRC) probability P with 10 neurons, (a) P = 0.0 (one-dimensional regular ring with periodic boundary condition); (b) P = 0.5; (c) P= 1.0. Here we 
should mention that the interactions between neighboring neurons are bidirectional (shown by bidirectional arrowed lines), while the long-range 
drivings (LRDs) are unidirectional (shown by unidirectional arrowed lines). Time delays are only considered in these unidirectional LRDs. 
doi:10.1371/journal.pone.0096415.g001 



parameter region. Fig. 3(c) exhibits a completely synchronous 



R- 



10 



(7) 



as an order parameter to measure the degree of synchronization 
and the related transitions induced by time delay and long-range 
connection in Newman-Watts SWNNs. 

Results 

Time Delay Induced Synchronization Transitions 

In this part, we firstly investigate time delay induced synchro- 
nization transitions in Newman-Watts SWNNs at certain LRC 
probability. Figure 2 displays the dependence of synchronization 
parameters R {\0 samples for each t, depicted by black dots) and 
R (the average of Rs for 10 samples, depicted by red dots) on time 
delay t at P=\.0. Four distinct parameter regions have been 
revealed by synchronization parameters as time delay is increased. 
When time delay is small (t <2.6), synchronization parameters are 
all close to zero. It indicates that the states of individual neurons 
are significantly different and the whole network oscillates 
asynchronously at all (domain I in Fig. 2, called as asynchronous 
region). It means that small time delay has no effect on 
synchronization in delayed Newman-Watts SWNNs. A typical 
asynchronous spatiotemporal pattern is shown in Fig. 3(a) for 
T = 1.0. In the white regions, the nodes fire, while in the black ones 
they are quiescent. Time passes from left to right. Most of neurons 
in the network oscillate asynchronously and irregular spatiotem- 
poral dynamics is observed. As t is in the narrow region of 
[2.8,3-0], some synchronization parameters increase abruptly. It 
indicates that synchronous performance of neuronal network 
improves remarkably in some samples. A weak synchronization 
state for t = 2.8 is revealed in Fig. 3(b). The excitatory fronts are 
more ordered both in time and space. Time delay induced 
synchronization transition has been detected in Newman- Watts 
SWNNs. And we call this narrow parameter region as the 
transition region (domain II in Fig. 2, indicated by grey rectangle). 

When time delay is moderate (3.2<t<5.4, domain III in 
Fig. 2), synchronization parameters R and R jump to unity 
simultaneously. It implies that moderate time delay in information 
transmission can induce complete synchronization in Newman- 
Watts SWNNs. Therefore, synchronous region is defined in this 
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Figure 2. Time delay induced synchronization transitions. 

Dependence of synchronization parameters R 00 samples for each i, 
depicted by black dots) and R (the average of Rs for 10 samples, 
depicted by red dots) on time delay t at P= 1.0. Four distinct parameter 
regions, i.e., asynchronous region (domain I for small t), transition 
region (domain II for narrow region of time delay i, indicated by grey 
rectangle), synchronous region (domain III for moderate t) and 
oscillatory region (domain IV for large t, indicated by blue rectangle) 
are revealed. 

doi:1 0.1 371/journal.pone.009641 5.g002 

spatiotemporal pattern for t = 4.0. All neurons in the network fire 
simultaneously and damp to their rest state together. As time delay 
is further increased (t>5.6), to our surprise, desynchronization 
occurs in Newman- Watts SWNNs. A distinct new parameter 
region, composed by asynchronous state, weak synchronization 
and complete synchronization, has been discovered. And oscillat- 
ing behaviour of the order parameter is detected. Accordingly, we 
call this parameter region as the oscillatory region (domain IV in 
Fig. 2, indicated by blue rectangle). Fig. 3(d) displays a typical 
desynchronized spatiotemporal dynamics in oscillatory region at 
T = 7.0. Large time delay can effectively improve synchronization 
in the beginning (can be indicated by the second excitatory front in 
Fig. 3(d)). However, the ordered excitatory front degenerates and 
desynchronization occurs as the system evolves. Finally, asynchro- 
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Figure 3. Space-time plots of u for different time delay t at P = 1.0. 

(a) T=1.0 (asynchronous state), (b) t = 2.8 (weak synchronization), (c) 
1 = 4.0 (complete synchronization), (d) t = 7.0 (desynchronized state). The 
figures are plotted in greyscale from black (lowest value at 0.0) to white 
(highest value at 1.0). And this greyscale will be used throughout this 
paper. 

doi:1 0.1 371/journal.pone.009641 5.g003 

nous state is obtained in oscillatory region. According to the results 
shown in Fig. 2, we can conclude that moderate time delay is 
needed for synchronization in delayed Newman-Watts SWNNs. 

For further investigating the synchronous oscillations, the 
dependence of oscillation period T on time delay t in synchronous 
region is shown in Fig. 4(a). It is seen that synchronization 
oscillation period is monotonously increased with time delay. And 
approximate linear relationship is revealed. However, a time 
difference between T and t can be detected. To explain the above 
phenomenon, time series u of neurons 79 (shown by black curve), 
78 and 80 (two neighboring neurons of 79, shown by green and 
yellow curves) and 42 (the LRD neuron of 79, shown by red curve) 
of Fig. 3(c) are shown in Fig. 4(b). The blue dashed curve denotes 
time series u of neuron 42 with time delay translation. The pink 
line indicates excitation threshold. From Fig. 4(b) we can find that 
synchronization oscillation period T is composed by time delay T 
and excitation time Ie- That's why there exists a time difference 
between synchronization oscillation period and time delay. 

The mechanism of synchronous oscillations can also be 
explained by Fig. 4(b). As complete synchronization is achieved 
in delayed Newman- Watts SWNNs, all neurons can excite 
simultaneously and damp to their rest state together, oscillate just 
as a single cell (can be indicated by the overlap of the four solid 
curves). Since time delays exist in LRCs, neurons can be excited 
synchronously again by their corresponding delayed LRDs (can be 
indicated by the black solid and blue dashed curves). Synchronous 
oscillations can self-sustain in delayed Newman-Watts SWNNs in 
this manner (such as the two excitation periods shown in Fig. 4(b)). 
However, due to the existence of refractory period for excitable 
neuron, a minimal time delay Tmin is needed for LRDs sustaining 
synchronous oscillations. Accordingly, complete synchronization 
can emerge in delayed Newman-Watts SWNNs as t > Tmin ■ Based 



on the results shown in Fig. 2, we can fmd Tmin^^-^ under current 
parameter settings. Now the transition from non-synchronization 
to complete synchronization can be explained as follow: For small 
time delays (i.e., KTmw), LRDs can not occupy the whole 
network entirely and simultaneously due to the existence of 
refractory period for excitable dynamics. Neurons in the network 
are mostly driven by their neighbors. As a result, zigzag excitation 
fronts (i.e., asynchronous spatiotemporal patterns) are obtained. As 
Tfnin is reached, LRDs can dominate the neuronal network 
absolutely, and complete synchronization can emerge in delayed 
Newman-Watts SWNNs. 

From the above discussion, we can find that LRDs play a key 
role in synchronization transition. To qualitatively investigate the 
effects of LRDs on the spatiotemporal dynamics obtained in 
delayed Newman-Watts SWNNs, the LRD proportion p is used, 
which can be calculated as: 



where Nlrd is the total number of neurons driven by LRDs. The 
evolvement of LRD proportion p between adjacent intervals for 
different time delay t (corresponding to Figs. 3(a)-d)) is shown in 
Fig. 4(c). As time delay is small (t = 1 .0, below Tmin, shown by black 
squares), p increases slightly at first and then tends to 0.5. It 
indicates that the neuronal network is governed by local and long- 
range drivings together. Accordingly, irregular asynchronous 
spatiotemporal dynamics of Fig. 3(a) is obtained. When t is in 
the transition region (t = 2.8, close to T^in, shown by green 
triangles), LRD proportion p increases abruptly, but can never 
reach 1.0. It means that most of neurons in the network are 
sustained by LRDs, and can fire simultaneously. However, few 
neurons are still excited by their corresponding neighbors. 
Therefore, weak synchronization can be observed. As synchronous 
region is reached (t = 4.0, beyond Tmin, shown by red dots), LRD 
proportion p jumps to unity rapidly. With the help of moderate 
time delay, LRDs can suppress neighboring interactions to 
dominate the system entirely. All neurons in the network can be 
excited by their corresponding LRDs simultaneously, and 
complete synchronization can emerge in delayed Newman-Watts 
SWNNs. When time delay is large (t = 7.0, also beyond Tmin, 
shown by blue diamonds), p increases abruptly at first and goes 
through a peak, then deceases monotonously, and finally tends to 
0.5. It means that LRDs can take effect so long as Tmin is reached. 
And LRDs can dominate the neuronal network in the beginning 
and weak synchronization such as the second excitatory front of 
Fig. 3(d) can be achieved. However, LRD loses its predominance 
as the system evolves. It may be caused by the too long resting time 
which can increase the chance for neighboring interactions. As a 
result, the ordered excitatory front degenerates and desynchroni- 
zation occurs in the oscillatory region. So we consider that too 
large time delay may be harmful for synchronization to a certain 
degree. 

Based on the above discussion, we can infer that the mechanism 
behind spatiotemporal dynamics obtained in delayed Newman- 
Watts SWNNs is the competition between LRDs and neighboring 
interactions. This kind of competition is caused by inhomogeneity 
in information transmission between neighboring and long-range 
interactions of the present model. More importantly, the 
competition results, which will decide the spatiotemporal dynamics 
in the network, are largely dependent on time delays. Therefore, 
we can expect that the LRD proportion is also a good indicator to 
study the synchronization transitions in delayed Newman- Watts 
SWNNs. The dependence of LRD proportion p {10 samples for 
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Figure 4. Dynamical analysis of synchronous oscillations and time delay induced synchronization transitions, (a) Dependence of 
oscillation period T on time delay t in synchronous region, (b) Time series u of neurons 79 (shown by black curve), 78 and 80 (two neighboring 
neurons of 79, shown by green and yellow curves) and 42 (the LRD neuron of 79, shown by red curve) of Fig. 3(c). The blue dashed curve denotes 
time series u of neuron 42 with time delay translation. The pink line indicates excitation threshold. The oscillation period T is composed by time delay 
T and excitation time Ie- (c) The LRD proportion p between adjacent intervals for different time delay t (corresponding to Figs. 3(a)-(d)). (d) 
Dependence of LRD proportion p [^0 samples for each t, depicted by black dots) and p (the average of ps for 10 samples, depicted by red dots) on 
time delay t. The four distinct parameter regions can also be revealed by LRD proportion clearly. 
doi:1 0.1 371/journal.pone.009641 5.g004 



each T, depicted by black dots) and p (the average of Ps for 1 0 
samples, depicted by red dots) on time delay t is shown in Fig. 4(d). 
The four distinct parameter regions are revealed by LRD 
proportion clearly. Moreover, we can also find that moderate 
time delay can help LRDs to beat neighboring interactions to 
dominate the network absolutely. The conclusion that moderate 
time delay is needed for synchronization in delayed Newman- 
Watts SWNNs is further verified. 

LRC Induced Synchronization Transitions 

From the above understanding we can find that LRDs play an 
important role in deciding the spatiotemporal dynamics. There- 
fore, a detailed study on LRC induced synchronization transitions 
needs to be taken in delayed Newman-Watts SWNNs. Fig. 5(a) 
displays the dependence of synchronization parameter R on LRC 
probability P for different time delay T. For small time delay 
(t=1.0, below Tmin^ shown by black triangles), LRDs can not 
occupy the system due to the existence of refractory period. As a 
result, LRCs have no effect on synchronization transitions in 
asynchronous region. When time delay is in transition region 
(t = 2.8, close to Xmin^ shown by pink squares), few LRDs can 
occupy the neuronal network under this circumstance. Therefore, 
lots of LRCs are needed to slightly improve the synchronization. 
For moderate time delay (t = 4.0, beyond Xmin^ shown by red dots). 



LRDs can suppress neighboring interactions to dominate the 
system entirely. Consequentially, synchronization in delayed 
Newman-Watts SWNNs can be enhanced remarkably by increas- 
ing LRC probability P. For large time delay (t = 7.0, also beyond 
Xmin^ shown by blue diamonds), synchronization of delayed 
Newman -Watts SWNN improves as LRC probability increases. 
However, as we have identified, too large time delay can increase 
the chance for neighboring interactions and is harmful for 
synchronization to a certain degree, oscillating behaviour of the 
order parameter can be observed. Fig. 5(b) displays the 
dependence of synchronization parameter R on time delay t for 
different LRC probability P. An optimal time delay interval is 
needed to enhance the synchronization for Newman-Watts 
SWNNs. The centers of optimal time delay interval are all around 
4.5 and are largely independent of LRC probability. The width of 
optimal time delay interval broadens as LRC probability increases. 

To give more intuitive understanding on LRC induced 
synchronization transitions in delayed Newman-Watts SWNNs, 
space-time plots of w for different LRC probability P at t = 4.0 is 
given in Fig. 6. Remarkable enhancement of synchronization 
induced by LRCs in delayed Newman-Watts SWNNs is revealed 
obviously. Besides the asynchronous state (P = 0.30 for Fig. 6(a)), 
weak synchronization (P = 0.70 for Fig. 6(b)) and complete 
synchronization (P=1.00 for Fig. 6(d)), another new synchroni- 
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Figure 5. LRC induced synchronization transitions, (a) Dependence of synchronization parameter R on LRC probability P for different time 
delay t. (b) Dependence of synchronization parameter R on time delay t for different LRC probability P. 
doi:1 0.1 371/journal.pone.009641 5.g005 



zation mode has been found at P = 0.96 and is shown in Fig. 6(c). 
From visual assessment, we guess this kind of new synchronization 
mode is the lag synchronization. To test our idea, the similarity 
function is introduced, which was proposed to detect lag 
synchronization [5]. It is numerically calculated as: 

[<V79(0'><V78(0' 

Here V79 and Vyg are the time series v of neurons 79 and 78 of 
Fig. 6(c). And is the time shift. Fig. 7(a) displays the dependence 
of similarity function S on time shift T^. The minimal value of S 




1 O 30 SC3 t3C3 



t (time units) 

Figure 6. Space-time plots of u for different LRC probability P 

at T = 4.0. (a) P = 0.30 (asynchronous state), (b) P = 0.70 (weak 
synchronization), (c) P = 0.96 (lag synchronization), (d) P=1.00 
(complete synchronization). 
doi:1 0.1 371 /journal.pone.009641 5.g006 



appears at t^q =0-363, which indicates the lag synchronization 
between neurons 79 and 78. Fig. 7(b) shows the projection of the 
attractor on the time shifted plane {vjs{t-\-Tso), ^79(0)- It 
demonstrates that the state of neuron 79 is delayed in time with 
respect to neuron 78. Accordingly, lag synchronization has been 
confirmed in delayed Newman-Watts SWNN. To explain the 
mechanism of lag synchronization, time series u of neurons 79 
(without LRC, shown by black curve), 78 and 80 (two neighboring 
neurons of 79, shown by red and blue curves) of Fig. 6(c) are 
shown in Fig. 7(c). And the red dotted and blue dashed curves 
denote time series u of neurons 65 and 93 (the two LRD neurons 
of 78 and 80) with time delay translation, respectively. As LRC 
probability P is a little less than 1.0, some neurons in network will 
have no LRCs due to finite connection probability. All neurons 
without LRCs must be driven by their neighbors. And these 
neighboring neurons are excited by their corresponding delayed 
LRDs. The successive driving relationship is revealed in Fig. 7(c). 
And lag synchronization between neurons without LRCs and their 
corresponding neighbors is identified. Therefore, we can observe 
lag synchronization in delayed Newman-Watts SWNNs as LRC 
probability P is a little less than 1.0 at moderate time delay. 
Fig. 7(d) exhibits the LRD proportion p between adjacent intervals 
for different LRC probability P at t = 4.0 (corresponding to 
Figs. 6(a)-(d)). Anticipated LRD proportions can be quickly 
approached so long as time delay is moderate. And large numbers 
of LRCs are needed to dominate the network for synchronization 
under this circumstance. 

According to the results obtained in this part, the conclusion 
that moderate time delay can help LRDs to dominate the network 
has been verified again. And large numbers of LRCs are needed 
for synchronization under this circumstance. Therefore, the two 
necessary conditions, moderate time delay and large numbers of 
LRCs, are exposed explicitly for synchronization in delayed 
Newman-Watts SWNNs. 

The Combined Effects on Synchronization Transitions 

To have a overall inspection of time delay and LRC induced 
synchronization transitions in Newman-Watts SWNNs, the 
contour plot of synchronization parameter R in the plane {t,P) 
is revealed in Fig 8. The color intensity denotes the synchroni- 
zation degree in delayed Newman- Watts SWNNs. Specifically, 
lighter color representing larger synchronization parameter, which 
indicates higher degree of synchronization. The four distinct 
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Figure 7. Dynamical analysis of lag synchronization and LRC induced synchronization transitions, (a) Dependence of similarity function 
S on time shift t^. The minimal value of S appears at t^o = 0.363, which indicates the lag synchronization between neurons 79 and 78 of Fig. 6(c). (b) 
Projection of the attractor on the time shifted plane {vjsit + Tso), v^git)). It demonstrates that the state of neuron 79 is delayed in time with respect to 
neuron 78. (c) Time series u of neurons 79 (without LRC, shown by black curve), 78 and 80 (two neighboring neurons of 79, shown by red and blue 
curves). The red dotted and blue dashed curves denote time series u of neurons 65 and 93 (the two LRD neurons of 78 and 80) with time delay 
translation, respectively. Lag synchronization is discovered in delayed Newman-Watts SWNN and the mechanism is also revealed, (d) The LRD 
proportion p between adjacent intervals for different LRC probability P (corresponding to Figs 6(a)-(d)). 
doi:1 0.1 371/journal.pone.009641 5.g007 



parameter regions, i.e., asynchronous region, transition region, 
synchronous region and oscillatory region at certain LRC 
probability P=\.0 are exposed clearly. And the remarkable 
enhancement of synchronization transitions induced by LRCs 
under moderate time delay is also indicated explicitly. From Fig 8 
the optimal combinations of time delay and LRC probability on 
synchronization transitions in delayed Newman-Watts SWNNs 
are revealed intuitively, which may has a useful impact for actual 
biological systems. 

The Universality of Time Delay Induced Synchronization 
Transitions 

In order to test the universality of time delay induced 
synchronization transitions, heterogeneous Newman-Watts 
SWNNs are considered. Diversity is introduced to system 
parameter b, i.e., the values of bi are different in the network. 
And it satisfies the following Gaussian distribution: 



SWNNs. Fig. 9(a) shows the dependence of synchronization 
parameter R on time delay t for different diversity a at LRC 
probability P=\.0. Although synchronization transition becomes 



<bi> =b,< (bi - 



-b)>=<j'd(i-j). 



(10) 




The value of ^ is fixed at 0.07. Here cr is the standard deviation 
of the Gaussian probability distribution of system parameter b. It 
indicates the strength of the diversity in delayed Newman- Watts 



Figure 8. The combined effects on synchronization transitions. 

Dependence of synchronization parameter R on time delay t and LRC 
probability P. 

doi:1 0.1 371/journal.pone.009641 5.g008 
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Figure 9. The universality of time delay induced synchronization transitions, (a) Dependence of synchronization parameter R on time 
delay t for different diversity a. (b) Dependence of synchronization parameter R on time delay t for the new coupling. 
doi:1 0.1 371 /journal. pone.009641 5.g009 



less profound as diversity increases, similar time delay induced 
synchronization transitions can be observed in heterogeneous 
Newman-Watts SWNNs. More importantly, all these synchroni- 
zation transitions appear at approximatively same t. It indicates 
that time delay plays a significant role in synchronization 
transitions in Newman-Watts SWNNs. Moreover, to further test 
the generality of our findings, the following new coupling form is 
used: 



dUi(t) 1 nr /X ^/(0 + ^ 

— ^ = - -Ui(t)[Ui(t)-\][Ui(t)- 

at e a 



(11) 



Wu + Ck 



-Ui{t)), 



where 



- [ui- 1 (0 + W/+ 1 (0 - 2ui{t)\ + [uj{t - t) - Uiit)] . (12) 



This type of coupling has been widely used in neural models and 
excitable complex networks. In simulations, we set Cj^ = 0.5. 
Fig. 9(b) displays the dependence of synchronization parameter R 
on time delay t. Similar time delay induced synchronization 
transitions can also be observed for the new coupling form. Now 
we can conclude that time delay induced synchronization 
transitions in Newman-Watts SWNNs are a robust phenomenon. 
The results revealed in present paper are universal. 

Conclusions 

In conclusion, time delay and long-range connection induced 
synchronization transitions in Newman-Watts small-world neuro- 
nal networks are systematically investigated by synchronization 
parameter and space-time plots. We have found four distinct 
parameter regions, i.e., asynchronous region, transition region. 



synchronous region and oscillatory region, at certain LRC 
probability P=1.0 as time delay is increased. Interestingly, 
desynchronization and oscillating behaviour of the order param- 
eter are observed in oscillatory region. More importantly, the 
mechanisms of synchronous oscillations and the transition from 
non-synchronization to complete synchronization are discussed. 
We consider the spatiotemporal patterns obtained in delayed 
Newman-Watts SWNNs are the competition results between long- 
range drivings and neighboring interactions. And our point of view 
has been verified by LRD proportion, which can also reveal the 
four distinct parameter regions clearly. In addition, for moderate 
time delay, the synchronization of neuronal network can be 
enhanced remarkably by increasing LRC probability. Further- 
more, lag synchronization has been found between weak 
synchronization and complete synchronization as LRC probability 
P is a little less than 1.0. Finally, the two necessary conditions, 
moderate time delay and large numbers of LRCs, are exposed 
explicitly for synchronization in delayed Newman-Watts SWNNs. 

As we know that synchronization transitions in neuronal 
networks are very important issues in related research fields and 
are associated with some specific physiological functions. A 
systematical investigation of synchronization transitions induced 
by time delay and long-range connection is expected to be useful 
both for theoretical understandings and practical applications. The 
results obtained in the present paper are universal. Similar time 
delay induced synchronization transitions can also be observed for 
heterogeneous Newman-Watts SWNNs and the new coupling 
form. We do hope that our work will be a useful supplement to the 
previous contributions and will have a useful impact in related 
fields. 
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